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In this work we study a disordered binary Ising model on the square lattice. The model system 
consists of two different particles with spin- 1/2 and spin- 1, which are randomly distributed on 
the lattice. It has been considered only spin nearest-neighbor exchange interactions with J > 0. 
This system can represent a disordered magnetic binary alloy A x Bi- x , obtained from the high 
temperature quenching of a liquid mixture. The results were obtained by the use of Monte Carlo 
simulations for several lattice sizes L, temperature T and concentration x of ions A with spin- 
1/2. We found its critical temperature, through the reduced fourth-order Binder cumulant for the 
several values of the concentration x of particles (spin-1/2, spin-1), and also the magnetization, the 
susceptibility and the specific heat as a function of temperature T. 



I. INTRODUCTION 

Magnetic properties of binary site-substitutionally disordered Ising models have been receiving considerable atten- 
tion from both theoretical and experimental point of view (If. In these systems, two different types of magnetic ions 
(denoted by A and B) are randomly distributed on a lattice representing a magnetic binary alloy A X B±_ X , which is 
suddenly frozen from high (liquid state) to low temperatures (solid state). Moreover, this systems may exhibit a rich 
diagram phase. 

Much of the theoretical work has assumed the two magnetic ions having the same spin-1/2 value. These models have 
been investigated by mean- field approches [H-Q and Monte Carlo simulation Q. On the other hand, less attention has 
been given to binary random-site systems where the constituints have different spin values. Thefore, it is interesting 
to investigate a system as the binary random-site Ising model with one of the constituints having spin-1 and the other 
having spin-1/2. These systems were investigated by mean- field theories [6H8| and Monte Carlo simulation [9j. In 
the reference Q the authors studied the system with half of the lattice with spin-1/2 and another half with spin-1 
randomly distributed. 

In this work, we used the mixed-spin Ising model approach for a binary alloy. The results were obtained by the use 
of Monte Carlo simulations for several lattice sizes L, temperature T and for several values of the concentration x of 
ions A with spin-1/2. We found its critical temperature, through the reduced fourth-order Binder cumulant for the 
several values of the concentration x of the particles (spin-1/2, spin-1), and also the magnetization, the susceptibility 
and the specific heat as a function of temperature T. 

The paper is organized as follows: in Section II, we describe the disordered binary Ising model and we define some 
observables of interest. In Section III, we present some details concerning the simulation procedures and the results 
obtained. Finally, in the last Section, we present our conclusions. 

II. MODEL AND SOME OBSERVABLES OF INTEREST 

The configurational energy of disordered binary Ising model may be described by the following Hamiltonian 

TL = - J ^{<y ,a,rr, + (1 - Ci)(l - Cj)SiSj + Ci(l - CjfaSj + cj(l - c l )S l a J }, (1) 

(id) 

where the spin variables assume the values Si = ±1,0 and <jj = ±1/2, and the nearest-neighbor interaction is 
ferromagnetic, J > 0. We associated to each site i of the lattice a occupation variable Cj, so that Cj = 1, if the site 
is occupied by a particle with spin a = 1/2 and Cj = 0, if it is occupied by a particle with spin S = 1. The sites are 
occupied independently by the two particles, with a probability distribution defined as 
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P(ci) = [(1 - x)8 Cjt0 + x8 Cjt 



i=l,2,...,JV 



(2) 



where x is the concentration of spin-1/2 and (1 — x) is the concentration of sphi-1. 
We calculate the following thermodynamic quantities per site: the magnetization mi 



the energy El 



1 

TV 



JV 



{(i - c^Si + ciOi} 



(3) 



Er 



the susceptibility \l 



and the specific heat Cl 



N 



XL=N([(m 2 L )] - [<m L ) 2 ]) 



c L = N ([{El)] [(El) 2 }) . 
To find the critical point, we used the fourth-order Binder cumulant (loL [Tlj. 

[«)] 



U L = 1- 



(4) 



(5) 



(6) 



(7) 



In the expressions above [■ • •] denotes average over the samples of the system, and (• • •) denotes the thermal average. 
N = L 2 is the total number of particles. 



III. PROCEDURES AND RESULTS 



In order to determine the critical behavior of this disordered binary Ising model we employed Monte Carlo simulation 
techniques. We consider a square lattice of linear size L, with values of L ranging from L — 12 to L = 48, and we 
applied periodic boundary conditions. 
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Figure 1: Snapshots of the spin configurations are given for two different concentrations x of spin- 1/2 and two different values 
of temperatures T. For x = 0.25, T = 0.5 (a) and T = 1.5 (b); x = 0.75, T = 0.5 (c) and T = 1.5 (d). spin-1/2 (spin-1) are 
represented by yellow (black). The temperature is measured in units of J/fcs. 

We prepared the system with the spins randomly distributed on the lattice. The concentration x of spin-1/2 and 
(1 — x) spin-1, is fixed. However, each spin-1/2 can have as its nearest-neighbor spins of the type spin-1/2 or spin-1 
and vice versa. Each trial change of a spin state on the lattice is accepted according to the Metropolis prescription 
(l2^ . oj(a —> a') = min[l, exp(— (3AE)], where AE is the local energy change [E a i — E a ) resulting from changing 
the state of a random selected spin from a to a' state, and j3 = l/fc^T. To reach the equilibrium state we take for 
guarantee at least 1 x 10 6 MCs (Monte Carlo steps) for all the lattice sizes we studied. Then, we take more 5 x 10 5 
MCs to estimate the average values of the quantities of interest. Here, 1 MCs means L 2 trials to change the state 
of a spin of the lattice. The average over the disorder was done by using 100 independent samples for lattices in the 
range 12 < L < 48. 

In Fig. 1, we have a snapshot of the concentrations of spin-1/2 and spin-1, and they were made for two different 
concentrations x — 0.25 and x = 0.75 of spin-1/2. In both cases, it has been analyzed the different concentrations at 
low temperatures (T < T c , T = 0.5) as shown in Fig. 1(a) and Fig. 1(c), and at high temperatures (T > T c , T = 1.5) 
as shown in Fig. 1(b) and Fig. 1(d). The temperature is measured in units of J/ks throughout the paper. 
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Figure 2: Thermodynamic quantities per site as a function of temperature T, and for several values of the concentration x 
of spin-1/2. (a) Magnetization m, (b) susceptibility x, ( c ) specific heat C and (d) fourth-order Binder cumulant U. The 
concentration x of spin-1/2 is x — 1.0 (gray line), x = 0.75 (blue line), x = 0.50 (green line) , x = 0.25 (red line) and x = 
(black line). The temperature is measured in units of J/ks- 

All the results of the simulations presented here were realized for a concentration x of spin-1/2 fixed at x = 0, 
0.25, 0.50, 0.75, 1.0. In Fig. 2, we displayed the behavior of the magnetization (Fig. 2(a)), susceptibility (Fig. 2(b)), 
specific heat (Fig. 2(c)) and fourth-order Binder cumulant (Fig. 2(d)) as a function of temperature T, and several 
values of the concentration x of spin-1/2. For x = 0, we have the critical behavior of the Blume-Capel model with 
the term of anisotropy equal zero (D = 0) fl5l. [l6|. Nevertheless, for x = 1.0, we have the case of the pure Ising 
model. The critical behavior for x = 0.50, i. e., half of the lattice with spin-1/2 and another half of the lattice with 
spin-1 randomly distributed were studied in the reference We can observe that the magnetization vanishes with 
the increase of temperature T. When the concentration x of spin-1/2 increases (x — > 1) the magnetization vanishes at 
different values of critical temperatures T c . The critical temperature decreases with the increase of the concentration 
x of spin-1/2. These can be clearly observed by the shifting in the susceptibility and specific heat peak (see Fig. 2(b) 
and 2(c)). 
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Figure 3: Magnetization mz, as a function of temperature T for various lattice sizes L as indicated in the figure. The temperature 
is measured in units of J/fcg. 

The finite-size effects of the magnetization were studied for the concentration x of spin-1/2 fixed at x = 0.75 and 
for various system sizes L. In Fig. 3, it is shown the behavior of the magnetization as a function of temperature T 
for several system sizes L. We can observe in this figure, that the magnetization tul vanishes with the increase of 
temperature T, indicating the existence of a phase transition. 
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Figure 4: Fourth-order cumulants Ul as a function of temperature T for various lattice sizes L as indicated in the figure. The 
temperature is measured in units of J/Ub- 

So as to study the phase transition in more details, we used the fourth-order Binder cumulants Ul intersection 
method to determine the value of temperature at which the transition occurs. According to the theory of finite-size 
scaling for continuous phase transitions, the finite-size behavior is governed by the ratio L/£, where £ is the correlation 
length. The scaling relation for the fourth-order cumulant shows that, at critical temperature, where the correlation 
length is infinite, all the curves must intercept themselves at a single point, since L/£ is zero for all the sizes L. To 
find the critical temperature, we displayed in Fig. 4 the cumulants Ul(T) versus temperature T, for several system 
sizes L and for x = 0.75. Our estimate for the dimensionless critical temperature is T c = (0.71 ± 0.01) J/ks- 

The susceptibility \l as a function of temperature T is shown in Fig. 5. For finite systems \l presents a peak 
around the critical temperature T c , which grows in height with the increase of the system size. 
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Figure 5: Susceptibility xl asa function of temperature T for various lattice sizes L as indicated in the figure. The temperature 
is measured in units of J/fcg. 

In Fig. 6, we also present the measurements of the specific heat. The peak observed in the curves of the specific 
heat exhibit a weak system size dependence compared to the susceptibility peak. The position of the specific heat 
and of the susceptibility peaks can be defined at a pseudocritical temperature T max (L). The T max (L) approaches T c 
when L — > oo [13| . 
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Figure 6: Specific heat cl as a function of temperature T for various lattice sizes L as indicated in the figure. The temperature 
is measured in units of J/fcg. 

We have also calculated the critical temperature T c for several different values of the concentration x of spin-1/2, 
as it can be seen in Fig. 7. For x = (only spin-1), we found T c = (1.69 ± 0.01) J/fc_e, that is the critical temperature 
of the Blume-Capel model with D = 0.0. Yet, for x = 1.0 (only spin-1/2), T c = (0.565 ± 0.003)4J/fc B which is the 
critical temperature of the Ising model. When x = 0.50 (half of spin-1/2 and half of spin-1) the critical temperature 
is T c = (1.00 ± 0.003) J/ks , therefore, the result is not according to the one of the reference^. The general features 
of the phase diagram show that the critical temperature change when concentration varies < x < 1.0. 
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Figure 7: Critical temperatures T c as a function of the concentration x of spin-1/2. The critical temperatures are measured in 
units of J/fcs. 



IV. CONCLUSIONS 

We have studied a disordered ferromagnetic binary Ising model on a square lattice. The model consists of two 
different types of particles with spin-1/2 and spin-1. These particles are randomly distributed on the square lattice, 
and we considered only nearest-neighbor interactions. We also calculated the critical temperature T c for several 
different values of the concentration x of spin-1/2. For x = (only spin-1), we found the critical temperature of the 
Blume-Capel model with D = 0.0. On the other hand, for x = 1.0 (only spin-1/2) we found the critical temperature 
of the Ising model. We have showed that the phase diagram present different critical temperatures for several values 
of the concentration x of spin-1/2. 
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